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Abstract 

The recent study by Bowers and Rajagopal on the possible struc- 
tures for the QCD crystalline color superconducting phase favors a 
cubic structure and allows for a first discussion of the QCD dynam- 
ics in the crystalline phase. The cubic structure leads to complete 
breaking of traslation invariance and to three phonon modes. We dis- 
cuss the expected form of the effective Lagrangian and the phonon 
dispersion laws assuming the proposed crystal structure. Besides the 
interest for QCD and the related possible relevance to compact stars 
these considerations may become useful to atomic systems exhibiting 
crystalline superfluidity. 



1 Introduction 

In a recent paper Bowers and Rajagopal have presented a comprehensive 
study of the different possible crystal structures for the crystalline color su- 
perconducting phase of QCD and concluded in favor of a face centered cubic 
condensate pattern. This study now allows discussion of the dynamics of 
QCD in the crystalline phase and this note is a first contribution in that 
direction. 
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The subject of the color superconducting phases of QCD [[| has only 
recently had its most extensive developments [|J. At asymptotically large 
values of the chemical potential for quark number /i, cold QCD with three 
degenerate flavors lives in the color-flavor locked (CFL) phase [j|] correspond- 
ing to condensation of quark pairs of opposite momenta and lying near the 
Fermi surfaces (BCS pairing). 

For finite strange quark mass the CFL phase is expected to persist within 
some range below asymptotic densities. Decreasing the value of /i, the ex- 
pected possibility becomes pairing between quarks of non vanishing total 
momentum, each quark momentum lying near the respective (now differing) 



Fermi surfaces @ [| [|T(J . Closeness to the Fermi surface will count for 
a lower free energy. Quantitative studies so far are mostly for two flavors and 
a difference in their chemical potential 5fi is usually introduced to induce the 
difference in the Fermi momenta. 

Similar possibilities of fermion pairing of non vanishing total momentum 
were already encountered in condensed matter (ferromagnetic systems with 
paramagnetic impurities) ||11||, so that the name LOFF phase is often used, 



from the initials of the authors of ref. 11 



In a LOFF phase translational and rotational invariance are sponta- 
neously broken. The most relevant question for crystalline superconductivity 
is that of the space dependence of the Cooper pair condensate. Preliminary 
studies have been based on a plane wave ansatz: a single plane wave with 
the finite momentum of the Cooper pairs, assuming that all the momenta 
spontaneously choose a unique direction. These studies (see ref || and UTTf ) 
have lead to a picture of the LOFF phase existing within a certain window 
of 5fi which starts with a first order phase transition and ends at the upper 
end with a second order transition. 

The authors of M develop a Ginzburg-Landau expansion near the second 
order transition. They introduce a multiple plane-wave ansatz and find that 
the favorite crystal structure is that of a face-centered cubic crystal. 

For the simple single plane wave structure, so far adopted as a provisional 
working hypothesis, we had constructed || the effective low energy phonon 
lagrangian and derived the phonon dispersion law. In that case only one 
phonon was expected from the partial breaking of translation invariance. 
For the cubic structure now derived in [0 one will have three phonons, due 
to the complete breaking of translation invariance. 

The work Jl| uses a simplified model with pointlike interaction and two 
massles quarks and, as we have said, a Ginzburg-Landau approximation. The 
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lowest value of the Ginzburg-Landau free energy corresponds to a crystal 
structure built up from eight plane waves directed to the corners of a cube. 
In spite of the approximations used one has the qualitative feeling that their 
proposed crystal structure would remain the favorite one also in a more 
general QCD setting (although not perhaps at asymptotic densities). 

In this note we shall discuss the expected form of the effective Lagrangian 
and phonon dispersion laws assuming the cubic crystal structure of 0. 

We have already mentioned the theoretical interest of color supeconduc- 
tivity in the crystalline phase as one of the general components of the QCD 
phase diagram. Physically, crystalline color superconductivity might be re- 
alized within compact stellar systems. In fact the quark chemical potential 
is expected to increase in penetrating towards the center of a compact star 
and other parameters, such as the strange quark mass, are also expected to 
evolve in such a way that one may reach the crystalline phase inside a core 
region of the star. In a rotating compact star the ensuing vortices may give 
rise to observable glitches, as discussed in ref. ||. 

The considerations developed in connection to crystalline superconduc- 
tivity in QCD, including those presented in this note, may also be useful in 
the field of atomic systems, in particular, as mentioned in ref. [ED], for gases 
of fermionic atoms below their Fermi temperature [ 12] [13|]. These systems 



present over QCD the advantage of a better parametric flexibility WA 



2 The LOFF phase of QCD 

According to |I| the favored crystalline structure in the LOFF phase of QCD 
should be a face-centered-cubic one. This is characterized by eight qj vectors 
with equal length q = n/a and directions corresponding to the eight vertices 
of the cube in Fig. 1: 

8 f 2iri } 
A(f) = A exp{2z(fj • x} = A ^ exp i (e 1 x 1 + e 2 x 2 + e 3 x 3 ) \ . (1) 

j=i n=± L a J 



3 




Fig. 1 - The figure shows some of the symmetry axes of the cube, denoted 
as C 2 , C 3 and C 4 (see the text). We show also the parameterization used for 
the coordinates of the vertices. 

In Fig. 1 some of the symmetry axes of the cube are shown: they are 
denoted as C4 (the three 4-fold axes), C3 (the four 3- fold axes) and C 2 (the 
six 2- fold axes). For more details on the symmetry group of the cube see the 



book by Hamermesh [15| . 

We observe that the coordinates at the vertices are nothing but the di- 
rections of the ijj vectors. The condensate (ttl) can be also written as follows: 



. . ~ r 2tx 2tt . 

A{x) = 2 A cos — (xi + x 2 + X3) + cos — [x\ — x 2 + £3) 
L a a 

2tt. . 2tt, , 

+ COS [Xi + X2 — X3) + COS l—Xx + X2 + X3) 

a a 



(2) 



This condensate breaks both translations and rotations. It is however in- 
variant under the discrete group Oh, the symmetry group of the cube. This 
can be seen by noticing that the previous expression is invariant under the 
following coordinate transformations 



R x : 


X\ 


R 2 : 


X\ 


R3: 


Xl 


I : 





Xl, X 2 — > X3, X3 — ) 
— X3, X 2 — ► X 2 , X3 
X 2 , X 2 -> -Xi, x 3 



-Xi, x 2 



-x 2 , 



x 2 , 

Xi, 

X3, 
X3 -> -X3, 



(3) 
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that is rotations of tt/2 around the coordinate axes, and the inversion with re- 
spect to the origin. Since one can easily check that the group Oh is generated 
by the previous 4 elements the invariance follows at once. 
We can then introduce 3 scalar fields $W such that 

9-7T 

($ W )o = =^Xi (4) 

These fields should be defined on the lattice sites, but for the purpose of 
writing the effective action for such fields, we are interested only in the low 
energy behavior p <C 1/a. Therefore we will treat them as continuous fields. 
The coupling of the Fermi fields generated by the condensate will be written 

as 

A^ T C^ ]T exp {z( ei $« + e 2 $ (2) + e 3 $ (3) )} (5) 

making the theory invariant under translations and rotations. These invari- 
ances are broken spontaneously in the vacuum defined by eq. (U). 

The effective lagrangian for the fields $^ has to enjoy the following sym- 
metries: rotational and translational invariance; Oh symmetry on the fields 
$W. The latter requirement follows from the invariance of the coupling fl5|) 
under the group Oh acting upon 

We then define 

<f>®(x) = —x i + <f>®( x )/f, (6) 
a 

with 

(0«(*)> O = O. (7) 

and / the decay constant of the phonon with dimension of an energy. It 
follows that, after the breaking, the phonon fields (jfo>(x) and the coordinates 
x l must transform under the diagonal discrete group obtained from the direct 
product of the rotation group acting over the coordinates and the Oh group 
acting over & % >(x). This is indeed the symmetry left after the breaking of 
translational and rotational invariance. 

In order to build up the effective lagrangian, we start noticing that if Xj 
are quantities transforming as the same representation of the $w's under Oh 
(a 3-dimensional representation), then all the invariant expressions can be 
obtained by the following three basic invariant expressions 

h{Xi) = Xf+Xf+Xf , J 4 (Xj) = X^X|+X|X|+X|X 2 , ie(Xi) = X\X\X\. 

(8) 
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The invariance is easily checked by noticing that these expressions are in- 
variant under the three elementary rotations Ri and the inversion I (see eq. 
(0)). Quantities transforming as 3>W are 



dt 



(9) 



We proceed now as in ref. ||, by noticing that the most general low-energy 
effective lagrangian, invariant under rotations, translations and Oh is given 
by 

c = t E ($ (i) ) 2 + £s(/ 2 (V$ (i) ), 7 4 (V$«), J 6 (V$W)) . (f 0) 

^ i=l,2,3 

The reason for allowing all the possible spatial gradients in this expansion 
is that the vacuum expectation value of is given by Xi and it is not in 

general small. 

Define the auxiliary quantities 



XAe) = — e W + eV0 
a 



(0 



such that 

Xi{l/f) = V$« , (12) 

where is a unit vector along the axis "i". The e parameter has been 
introduced to keep track of the spatial derivatives. Since we are interested at 
the second order in these derivatives we need the quantity (first order terms 
are total derivatives) 



1 d 2 C(e) 



2 de 2 



e=0 



E 

ij'=2,4,6 



d 2 C dhdL _ dCd 2 h 

+ E 



dlidlj de dt 



i=2,4,6 



dh dt 2 



(13) 



e=0 



where C(e) has been obtained by the substitution — > Xj(e). Since the 
derivatives of £ with respect to the invariant terms /«, evaluated at e = are 
just constant, using 



dt 



o = 2 E 



= 4 E <^ w , ^ 

£ =° >=W & 



£=0 



2 J] 



i=l,2,3 



(14) 
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and 



i d 2 i 2 

2 de 2 
ld 2 L 



2 de 2 
ld 2 h 



2 <9e 2 



e=0 



e=0 



e=0 



E 



(i)|2 



i=l,2,3 



4 E 0^8^ + 2 E |V0 W | 2 

i<j'=l,2,3 



i=l,2,3 



4 e + E |v</> w r. 

i<j'=l,2,3 i=l,2,3 



we see that the expression (|T3| ) depends only on the quantities 

E |V0«| 2 , E W (0 ) a . E . 

i=l,2,3 i=l,2,3 i<j'=l,2,3 



(15) 



(16) 



In conclusion one gets 

£ = ^ E E lw (i) i 2 -| E W w ) a -c E d^%^. 

1=1,2,3 i=l,2,3 1=1,2,3 i<i=l,2,3 

(17) 

This expression could have been guessed by noticing that after the breaking 
one still has an invariance under the group Oh acting at the same time upon 
the coordinates and the phonon fields (f>^\ 



3 Dispersion laws 

To analyze the phonon dispersion law one has to diagonalize the matrix 



/ a \p\ 2 + bp 2 cpip 2 



M 



cpxpz 



cpip 2 a\p\ 2 + bp 2 2 cp 2 p 3 

a \p\ 2 + bp\ 



1 2 i h^i 
\ Cp!p 3 Cp 2 p 3 

The dispersion law can be written as follows (p = np): 

E = v r (n)p 

where v 2 (n) are given by 



;is) 



(19) 



and u r are the three cubic roots defined by 



u r {n) = VA(n) + jB{ti) . 



(21) 



Here 



A(n) 
B(n) 



bJi — Jr 



/3 _ _ j 2 

■h 12 ^4 



b_ 



J4J6 



T 2 



(22) 



and 



J 4 



b 2 -& 



h(n) 



(b + 2c)(b-c) 2 h{n) 



(23) 



These solutions display the same symmetries of the effective lagrangian be- 
cause they depend on the directions n only through the elementary symmet- 
ric functions hin) = 1, hin), I^n). The corresponding eigenvectors can be 
computed in a straightforward manner; for example if 



rii 7^ , a — v 2 + bn\ 7^ , 
one gets the eigenvectors: 

/ 



v 2 + (b- c)n\ ^ 



(24) 



const, x 



V 



-n 3 c 



ni[a - v 2 + n\{b - c)] \ 

n 2 [a — v 2 . + n\(b — c)\ 

(a - v 2 )(nl + n 2 2 ) + 2(6 - c)n\n\ 
a — v 2 + bn\ / 



(25) 



The dispersion law for the phonons depends on three arbitrary parameters 
a, b, c. This follows from the symmetry of the cubic crystal and could be also 
inferred from the general theory of the elastic waves in crystals. As matter of 
fact, the phonon amplitudes (jP can be identified with the local deformations 
of the crystal with harmonic components 



(Pj = (f 0j e 



i(p-x-Et) 



(26) 
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The equation of motion for the deformation is 



P^T = ^ ( 27 ) 

where Oj\. is the stress tensor and p is the crystal density. The tensor djk is 
related to the strain tensor ue m : 

utm = ]z {de¥m + d m (fi + dnpid m ipi) , (28) 

by the relation 

Cjfc = ^iklmUim ■ (29) 

Kkim is the elastic modulus tensor and the equation for elastic waves can be 
written as 

P q^J = ^jkimdkdeipm , (30) 

from which the dispersion law follows: 

Det ( K pE 2 5 im - \iki m uimVkVi) = • (31) 

Owing to its symmetry properties, the elastic modulus tensor depends in 
general on 21 parameters; for the cube group Oh only three of them are 
independent [0, for example Aim, A1122, A1212, The identification with the 



parameters of the effective lagrangian is as follows: 

A1212 , Aim — A1212 A1122 + A1212 /or> N 
a = , = , c = . (62) 

P P P 



The group velocity 



1 (33) 



is a function of the direction of the wave propagation n but does not depend 
of the energy. For each fixed n there are three velocities that are in general 
all different. For comparison in an isotropic medium \v\ does not depend on n 
and two of the velocities (the transverse ones) coincide, while being different 
and smaller than the velocity of the longitudinal waves. 

Since the matrix M in ( |T8"D is real and symmetric, its eigenvalues are 
real as well. One has also to impose that v 2 r > 0; for this to happen all the 
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principal minors 

Ai = a + bnj , (z = 1, 2, 3) (34) 

A 2 = (a + bn^ia + bn 2 ) - c 2 n 2 n 2 , (i ^ j) (35) 

A 3 = a 2 (a + 6) + a(6 2 - c 2 )/ 4 (™) + (b + 2c)(b - c) 2 I 6 (n) (36) 

must be positive: 

A. 4 >0 (i = 1,2,3). (37) 

The constraint Ai > gives the conditions 

a > , a + 6 > . (38) 

The most stringent constraint A 2 > is obtained for nf = n 2 = 1/2, which 
gives 

b + 2a>\c\ . (39) 

Finally the most stringent constraint A3 > is obtained for I±(n) = 1/3 and 
Ie(n) = 1/27, i.e. when n is parallel to one of the 3- fold axes; since A3 is the 
product of the three eigenvalues, one gets 

b + 2c 

a + > . (40) 

3 

These constraints are summarized in Fig. 2. 

Let us specialize these results for momenta along the symmetry axes of 
the crystal i.e. along the three 4-fold axes, the four 3-fold axes and the six 
2-fold axes. 

We easily get for n along one of the three 4-fold axes, e.g. n = (1, 0, 0), 

v\ = a + b (41) 
v\ = v 2 3 = a (42) 

and the corresponding eigenvectors are 

6 = (43) 

while £3 are in the subspace spanned by (p^ and (p^ 3 \ The eigenvalues 
for the other two 3-fold axes with n = (0, 1, 0) and n = (0, 0, 1) are iden- 



tical to (|4"T| ) and ( |42| ) due to the cubic symmetry, while the corresponding 
eigenvectors are obtained by cyclic permutations. 
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Fig. 2 - The figure shows the regions where the square of all the velocities of 
propagation along the symmetry axes of the crystal are positive, in the plane 
(b, a) at fixed c. For c < only the area with horizontal lines is allowed. For 
c > both regions are permitted. 

For n along one of the four 3-fold axes, e.g. n = — /=(1, 1, 1) one has 

v3 

b + 2c 



3 

v 2 2 = v l = a + —^ (44) 



64(^ + ^ + ^) ( 45 ) 



with eigenvectors 

while £2, £3 are in the orthogonal subspace,i.e. the subspace spanned by 

- + x0 (2) + (1 - x)0 (3) , (46) 

with x a real parameter. 

The other three 3- fold axes are defined by the directions n = — —1,1), 

11 

n = — /=(— 1, 1) 1) an d ^ = T/^ - — ^ e Y can De obtained by the unit 

vector n = —^=(1,1,1) by applying transformations in O^, due to the cu- 
V3 

bic symmetry the eigenvalues are again given by (H) and the corresponding 
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eigenvectors are obtained by applying the same transformation to the eigen- 
vectors fl4"5|), (pESD ; for example for n = — t=(1, — 1, 1) one has 

v3 

6 = -L (0(D - 00) + 0(3)) (47) 

while £ 2) £3 are i n the orthogonal subspace spanned by 

(1) + x0 (2) + (1 - x)0 (3) . (48) 

Finally we consider the six 2-fold axes defined by the directions n = 
1 1 1 

— t=(1, ±1, 0), n = — 7=(1, 0, ±1) and n = — ^=(0, 0, ±1). By way of example 
v2 v2 v2 



we take 



The eigenvalues are 



n = -^(l,±l,0) . (49) 



o b + c , „ 

»i = (5°) 

9 6 — c . 

^2 = °+^- ( 51 ) 

^3 = a . (52) 



The eigenvector corresponding to the eigenvalue floTf ) can be obtained by 
®" Ci = ± (2) • (53) 



For the other two eigenvectors the conditions (24]) are not satisfied and one 
has to solve the eigenvalue equations by assuming n given by (fl9"D from the 
very beginning; one obtains 

6 = (3) ■ (54) 

The eigenvalues for n along the other 2-fold axes are identical, due to the 
cubic symmetry; exactly as in the preceding cases the corresponding eigen- 
vectors are obtained by (j53[) , (|54]) by applying the same transformations 
that change ( |4*9"D to the unit vectors of the other axes. All these results are 
summarized in Table 1 where we have introduced the unit vectors along the 
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coordinate axes, We observe again that in the case of propagation along 
axes of type C% there are three different propagation velocities, where for the 
other two cases the independent velocities are only two. The case of propa- 
gation along axes of type C3 is of particular interest since these directions are 
the same as the directions of the vectors q specifying the lattice. One of the 
two eigenvalues is doubly degenerate, and the non-degenerate eigenvectors is 
proportional to 

E Pi^- (55) 

8=1,2,3 

The simplest case is for a propagation along the coordinate axes (type C4), 
since in this case the matrix M is diagonal. Also in this case we have a 
doubly degenerate eigenvalue with eigenvectors the phonons corresponding 
to the plane perpendicular to the direction of propagation, whereas the non- 
degenerate eigenvector corresponds to the direction of propagation. It is 
interesting to notice that in all cases one of the eigenvectors is proportional 

to 



4 Conclusions 

Crystalline superconductivity arises from Cooper pairing between quarks of 
non vanishing total momentum, each quark momentum lying near the re- 
spective differing Fermi surfaces. Such a possibility (LOFF pairing) is known 
from study of certain condensed matter systems and recent theoretical stud- 
ies have made it as highly probable within the QCD phase diagram. LOFF 
pairing leads to spontaneous breaking of translational and rotational invari- 
ance. Depending on the particular crystal structure a number of phonons 
are the lowest energy manifestation of the breaking. 

A comparative study of the different possible crystal structures for the 
crystalline color superconducting phase of QCD has recently been published 
with the conclusion that a face centered cubic crystal is the favorite 
structure. Assuming this condensate pattern we have derived the expected 
form of the effective low energy Lagrangian and the dispersion law for the 
the set of the three phonons. The parameters a, b and c appearing in the 
effective lagrangian can be evaluated following the methods illustrated in ref. 



17fl . The results will be given elsewhere. 
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Momentum 


Eigenvalues 


Eigenvectors 






(a + 6) p 2 
ap 2 


0(0 

0(i), 0( fe ) (i^j ^ k) 






(3a + b - c) p 2 /3 
(3a + b + 2c) p 2 /3 


_0(1) +X 0(2) + (l_ x )0(3) 
0(1) + 0(2) + 0(3) 




(3a + b - c) p 2 /3 
(3a + b + 2c) p 2 /3 


0« + a:0( 2 ) + (1 - x)0 (3) 

_0(1) + 0(2) + 0(3) 


p = p(e (1) -e< 2 > + e< 3 >) 


(3a + b - c) p 2 /3 
(3a + b + 2c) p 2 /3 


0^ +^0^ - (l 

0(1) _ 0(2) + 0(3) 




(3a + 6 - c) p 2 /3 
(3a + b + 2c) p 2 /3 


0(1) + X 0(2) + (! + ^0(3) 
0(1) + 0(2) _ 0(3) 


c 2 


p — p (e ^ + e ^) 
i < j 


ap 2 

(2a + 6 - c) p 2 /2 
(2a + 6 + c) p 2 /2 


0(fc) 

_0« + 0(i) 

0« + 00) 


i < j 


ap 2 

(2a + b - c) p 2 /2 
(2a + 6 + c) p 2 /2 


0(fc) 

0(0 + 0(i) 
_0« + 0(i) 



Table 1 - Eigenvalues and eigenvectors for p along the symmetry axes. 
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